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Abstract 

A generalization of the classical secant condition for the stability of cascades of scalar 
linear systems is provided for passive systems. The key is the introduction of a quantity 
that combines gain and phase information for each system in the cascade. For linear one¬ 
dimensional systems, the known result is recovered exactly. 


1 Introduction 


An often-used tool in the analysis of biological feedback loops is the secant condition for linear 
stability; see the classical papers by Tyson and Othmer EH and Thron EDI, as well as the 
recent paper fSj. Consider a matrix of the following form: 


/ — a± 0 ... 0 — (3\ \ 

@2 — 0-2 • • • 0 0 

Vo 0 • • • P n -a n J 


where all ctj > 0 and all /3 ) > 0. Then, the secant condition states that the matrix is Hurwitz 
provided that: 


Pi • • • Pn 

OL\ . . . Oi n 



In essence, this says that a stable system with distinct real eigenvalues and no zeros tolerates 
negative feedback with a gain much larger than that provided by the small-gain theorem: the 
corresponding small-gain estimate would just have a “1” in the right-hand side. (The secant 
expression is always bigger than one. It is singular at n = 2 -which it should be, since then the 
matrix is always Hurwitz- and it equals 8 for n = 3, 4 for n = 4, and ~ 2.88 for n = 5, and 
tends monotonically to 1 as n —> oo. The bound is achieved exactly when all the aP s are the 
same.) The condition is useful because certain models of biological systems are stable for gains 
larger than those predicted by a simple application of the small-gain theorem. The secant takes 
advantage simultaneously of phase and gain information on the open-loop system. 


We provide here a generalization of the secant condition to cascades of output strictly passive 
(OSP) systems. We do so in such a manner that, in the special case in which each system is 
linear and one-dimensional, the classical result is recovered. (For simplicity, we restrict ourselves 
to systems with scalar inputs and outputs, but it is obvious how to generalize to more arbitrary 
systems, as long as inputs and outputs have the same dimension.) 
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The generalization is based on systematic use of a “gain” associated to OSP systems. It 
would seem that the use of this quantity might be useful for many other problems as well. 

This note is organized as follows. Section [21 introduces the basic concepts and states the 
main result, which is proved in Section ^ (the proof is actually very easy, given the definitions). 
Section [21 briefly mentions some extensions of the basic formalism, and Section [21 collects several 
facts concerning secant gains for the special case of linear systems. 


2 Notations, Definitions, and Statement of Main Result 

As usual, the extended space L 2 (0, oo) denotes the set of signals (thought of as time functions) 
w : [0, oo) — * R which have the property that each restriction wt = ^|[o,T] is in T 2 (0,T), for 
every T > 0. Given an element w E L 2 { 0, oo) and any fixed T > 0, one writes ||rc|| T for the L 2 
the norm of this restriction wt, and given two functions v,w E L 2 (0,oo) and any fixed T > 0, 
the inner product of vt and wt is denoted by ( v,w)t • In any Hilbert space, one defines the 
angle 9(v,w) E [0, 7 r] between two elements v,w by the formula 

Of ^ ( v ’ w ) 

COS 9(V,W) = |, |, -rr 

\\v\\ ||u;|| 

if v and w are nonzero, and zero otherwise. Given v, w E L 2 ( 0, oo) and any fixed T > 0, we will 
write 9t(v,w) instead of 9(vt,wt ), to denote the angle between the restrictions of the signals 
to [0, T]. 

We consider continuous-time finite-dimensional systems x = f(x,u), y = h(x) in the usual 
sense of control theory (e-g. B), with scalar valued inputs and outputs, and state space R ra , and 
assume always that the system is L 2 -well-posed, in the sense that for each u E L 2 (0, oo) and 
initial state x( 0 ) = 0 there is a unique solution x(-) defined for all t > 0 and the corresponding 
output y(t) = h{x(t )) is also in L 2 (0, oo). We call (u,y) is an input/output (i/o) pair of the 
system. 

We recall the standard notion of an output strictly passive (“OSP” for short) system, as 
given in textbooks such as Hues]. A system is OSP if there is some 7 > 0 such that, for 
every i/o pair ( u,y ), 

llz/ll t < 7 (u,y)r (1) 

for all T > 0. (Allowing an additive constant in the inequality is useful when dealing with 
arbitrary initial states. As we will study zero-state responses, we do not include a constant.) 

If a system is OSP, we call the smallest 7 as in Q the secant gain of the system, and denote 
it as 7 s . (There is a smallest such 7 , since the set of 7 ’s that satisfy Q is a closed set.) 

An equivalent definition of 7 $ is as the smallest 7 with the property that 

lll/ll T < TIMIt \\y\\T e T{u,y), 

or equivalently: 

\\y\\ T < 7 IMIt cos9 T (u,y) (2) 

for all T > 0 and all i/o pairs. Since m implies that {u,y)x > 0 for all i/o pairs and all T, 
for OSP systems we always think of the angle as lying in the interval [0, 7 r/ 2 ], and the cosine is 
nonnegative. 

The Cauchy-Schwartz inequality applied to 0 gives ||y|| T < 7 ||u|| T < 7 \\u\\ for all T > 0, 
so in particular y E L 2 if u E L 2 , and an OSP system necessarily has finite P 2 -induced (or 
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“H 0 o”) gain 7 ^ < (we remark later that this inequality is in general a strict one). Just 
as the L 2 gain is the supremum of the expressions ||?/|| T /||u|| T over all T and all i/o pairs 
with nonzero u, the secant gain is obtained by maximizing sec Ox(u,y) ||y|| T / ||u|| T , hence our 
terminology. 

If u E L 2 , so that also y E L 2 , taking limits in © gives 

II V\\ 2 < 7 (u,y). (3) 

Conversely, if u E L 2 => y E L 2 and © is true for all «£t 2 , then © holds. This is a routine 
exercise in causality, as follows. Pick any i/o pair (u,y) and any T > 0. Let v E L 2 be input 
which equals u on [0, T] and is zero for t > T, and z the output corresponding to v. Since 
v E L 2 , also 2 E L 2 . By causality, z restricted to [0, T] is the same as y restricted to [0,T], so 
(v,z) = ( u,y)x , and ||y|| T = ||z|| r . Therefore \\y\\^ = \\z\\^ < || z || 2 < 7 (v,z) = r y(u,y)T, and 
indeed © is verified. 

We wish to analyze the stability of the closed-loop system x = f(x,u — h(x)) obtained 
under negative unity feedback. Specifically, we study a cascade of n subsystems, as shown 
in the diagram in Figure Cl and subject to unity negative feedback. Such cascades appear 



Figure 1: Closed-loop system 


frequently in control theory as well as in biological applications, and, when components are 
one-dimensional, tend to have especially good dynamical properties such as the validity of the 
Poincare-Bendixson Theorem (|S3). We will assume that the z-th system has a secant gain 7 j, 
and we write y* for the output of the ith subsystem. We also assume well-posedness of the 
closed-loop. 

The main result is as follows: 

Theorem. Suppose that 

7172 ■■• 7 n < ( sec “) • 

Then the cascade is Testable: there is a number c so that 

11 Vn 11 T ^ C 11 U 11 rp 

for all input/output pairs in the cascade and all T > 0. 

Of course, this property implies as well that every | yj 11 T is bounded by some linear function 
of ||u|| T , and that the signals y* belong to L 2 if u E L 2 . 

For the special cases n = 1 and n = 2 (secant is infinite), we interpret the inequality in 
the theorem as saying that the condition holds for any possible values of the 7 j’s. For n = 2, 
therefore, the theorem is simply a restatement of the Passivity Theorem as given e.g. in na. 
Theorem 2.2.15, Part a (using only the input u). The Passivity Theorem also includes usually 
a statement (“Part b” in the citation) regarding the case in which the first system is OSP and 
the second one is only passive, meaning that only (u,y)x > 0 is known for all i/o pairs. We 
comment later on this fact. 

The assumption that the initial state of the cascade is x(0) = 0 is easy to dispose of, 
assuming appropriate reachability of the cascade, as routinely done in going from input/output 
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stability to state space stability, and Barbalat’s Lemma combined with either reachability or 
detectability arguments can be used to show convergence of internal states to zero. As an 
illustration, we state just one such corollary: 

Corollary. Suppose that the condition in the Theorem is verified, that the composite system 
shown in Figure n is zero-reachable and that each subsystem is input to state Testable. Then 
the system with no inputs (u = 0) has the property that all solutions converge to x = 0. 


3 Extensions 


We have formulated the results in terms of state-space systems only in order to be concrete. One 
could equally well consider arbitrary operators L 2 -> L 2 , or even just relations R on L 2 xL 2 , 
where an “i/o pair” is by definition any element of R, and define secant gain as the smallest 
number so that 0 holds for all T and all i/o pairs. Nor is it needed for the inputs and outputs 
to be scalar-valued; one may consider values on arbitrary Hilbert spaces, with inner product 
and norms taken pointwise in that space. More generally, functions of time are not required: 
one could consider an arbitrary Hilbert space H and simply ask that u and y belong to H. 
(To be precise, one needs a Hilbert space together with a resolution of the identity, in order 
to be able to be able to have a concept of “restriction” of u and y to subintervals; this is the 
formalism of resolution spaces developed in 0.) Even more generally, if one has a system in 
which inputs u and outputs y are known to lie in a specific subset S C H, then can be 
defined in terms only of i/o pairs that lie in S'; the validity of the main theorem is not affected, 
since it is just an algebraic statement about norms and inner products. 


Let us discuss a simple example of an operator defined only on subsets, which is of interest 
in biomolecular applications (“Michaelis-Menten kinetics”). Suppose that S is the set of all 
L 2 maps w : [0, oo) —► [—a, oo) with any fixed a > 0, and that we consider the function 
£ : [—a, oo) —> M given by 


£{r) 


Vr 

K + a + r 


( 4 ) 


(with K , V > 0 some constants) and the operator mh y defined on S, where y(t) = F(u)(t) = 
£(u(t)). This is an example of a “sector” nonlinearity. The analysis of sector nonlinearities is 
routine in passivity theory. The operator F is OSP and has 7 S = V/K, because we have, for 
all r G [—a, oo): 


W )] 2 = 


(since K + a + r > K), and thus 


12 
It 


V Vr 2 

K+a+rK+a+r 


< 


V Vr 2 V n , . 

<-= — r£(r) 

~ KK + a + r K v J 


f £{u{t)) 2 dt 

Jo 

V f T V 

k J 0 dt = 


so 7 S < V/K, and the equality is verified when u{t) = —a. 

Stability in the L 2 sense is only appropriate when dealing with equilibria associated to 
zero signals. However, the framework described here can be easily extended to more general 
situations. These extensions are of interest, particularly, when dealing with problems in biology 
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and chemistry, where quantities represent concentrations of substances, and hence are always 
nonnegative. We now describe briefly how this extension can be accomplished. 

Suppose that one wishes to study a system 

x = f(x,u ) 

V = h(x) 

under the feedback law u = — y, and that there is a steady state x* for this closed-loop system: 

f(x*,-h(x*)) = 0 

whose stability is of interest to analyze. We assume that the states x(t) evolve in some subset 
S of M ra , for example the positive orthant M” = {(aq,... ,x n ), aq > 0Vi}, and inputs u of the 
open-loop system take values on some set U. (In order for the closed-loop system to make 
sense, one should then have that —h(S) C U, of course.) We perform a change of variables 
z = x — x* and define the new system 

Z = g(z,v ) = f(z + x*,v — h{x*)) 
w = £{z) = h(z + x*) — h(x*), 

with states z(t) in the state-space {x — x*, x E S}, inputs v(t) in the input-value space {u + 
h(x*), u E 17}, and outputs w{t). Note that g(0,0) = 0. Applying the feedback v = —£(z) 
results in 

z = g{z, -£{z)) = f(z + x*,-h(z + x*)). 

Therefore, for each solution x(t) of x = f(x,—h(x)), the vector function z(t) = x(t) — x* 
satishes z = g(z,—£(z)), and conversely, each solution of the latter system arises from the 
former. Proving that solutions of x = f(x, —h(x)) converge to x* is then equivalent to proving 
that the solutions of the new system converge to the equilibrium z = 0. Thus we have reduced 
the analysis to the case treated in this paper. 

For example, suppose that we wish to study a positive system, that is, a system whose state 
state space is and inputs are also nonnegative. Furthermore, suppose that, as is often the 
case in biological feedback loops, one wishes to study an inhibitory feedback of the form 

M 

u = 77 - 

A T x n 

where M and K are some positive constants and x n is the nth coordinate of the state, that is 
to say, we have h(x) = —M/{K + x n ). In terms of the variables z, we have the output 

w = £(z) = h(z + x*) — h(x*) 

M M Vz n 

K + .T* K + (z n + x*) K -h x*+ z n 

which is the function in Q with a = x* and V = M/(K + x*). Since x n (t) is nonnegative, 
the state variable z n (t) takes values in [— x*, oo). Thus, we may view the closed-loop system as 
built from cascading the original system (which may itself be a cascade of several subsystems) 
with the static system “y = £(u)", which has % = V/K, and the previous analysis applies. 

This is all particularly simple for a linear system x = f(x, u ) = Ax + Bu. Positivity amounts 
to asking that all the off-diagonal entries of A as well as all entries of B are nonnegative 
(see e.g. in ei ). Since the system is linear and Ax* — Bh(x*) = 0, we have that g(z,v) = 
A(z + x*) + B(v — h(x*)) = Az + Bv , so the same open loop system results, except that now 
we are interested in the stability of z = 0. 
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4 Proof of Main Result 


Given an external input u, the solutions of the closed-loop system with initial state zero are so 
that the signals y t have the following properties: 

llz/ill t < n{u + yo,yi)T 
WvzWt ^ 72 {yi,y 2 )T 

hnWr < 7n (yn-l,y n )T 

for every T > 0, where we are writing y 0 = -y n . We expand (u + y 0 , yi) T = (u, yi) T + {yo, Vi)t, 
and use the Cauchy-Schwartz inequality for the first term, upper-bounding it by ||u|| j- \\yi || r . 
Replacing now each (yi-i,y t ) T by ||yj_i|| T ||yi|| T cos 9 T (yi-i, yi) and dividing by |;y ? ;|| T (assumed 
nonzero; otherwise, there will be nothing to prove), we have these estimates: 


2/1 T 

< 

7i ||2 /o|It cos <M2/o,2/i) + 7i INI t 

2/2 T 

< 

72 '2/1 t cos ^t(2/1)2/2> 

hn\\ T 

< 

7n \\yn-l\\ T cosS T(yn-l,yn) 


from which we conclude, by recursively substituting the estimates starting from the last one 
backward towards the first, that: 


112/«-11 x 1 — ^ ® II^IIt 


where 

a = 7172 • • • 7n c °s#r(2/i; 2 / 2 ) ■ ■ ■ ■ cos 0 T (y n -i,yn) 


and 


k = a cos d T (yoi Vi) ■ 


It is enough to show that k < 1, since then we can write (1 — /c) ||y n || r < a. ||it|| T , and therefore 
the result holds with c = a/{ 1 — k). Let us fix T and write 6{ := 0T(yi-i,yi) for i = 1,n. 
We must show, then, that 

/ 7T\ n 

cos 6\ ... cos 6 n < |^cos — J . (5) 


The angles 6i all lie in [0, vr/2], for each i = 2,..., n, since each system is OSP; thus cos 0- L > 0 for 
all such i. However, it is possible that cos 6\ < 0, since all that is known is that ( u+yo , yi)r > 0, 
not that ( yo,yi)T > 0. But if cos$i < 0, then (0) is true because the left-hand side is < 0 and 
the right-hand side is positive. So, in order to prove we may assume from now on that all 
9i e [0,vr/2]. 


We prove, more generally, this fact about Hilbert spaces: suppose given vectors vq,vi, ... ,v n 
such that (vi,Vi + i) > 0, and vq = —v n . Let 6i e [0, vr/2] be the angle between Vi-\ and u*. 
Then © holds. Intuitively, the property that the start and end vector are at angle n means 
that the consecutive vectors cannot be too close in angle, and therefore at least some of the 
angles must be large, and hence have small cosine, and the largest possible value is achieved 
when all angles are the same. 


To prove this general fact, without loss of generality, we may assume that all the are 
unit vectors (since only angles matter). Notice that Qi > n. This is because, for any three 
unit vectors, 6(u,v ) + 6(y,w) > 9(u,w), since we can view the angle as the geodesic distance 
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in a sphere, and apply the triangle inequality; inductively applied starting from vo, we get that 
— @( v o, v n) = 7T. Now, we have also this algebraic fact: 


cos 9\... cos 0 n < 



6\ + ... + 0 n 

n 


n 


which follows by noticing that the function f(x) = — In cos x is convex for x E [0, 7t/2) , applying 
Jensen’s inequality to obtain /QT0j/n) — (1 an d taking exponentials. Together 

with 0 % > 7 r, using that 7r jn < {0\ + ... + 0 n )/n < ir/2 < 7r (recall that each 0i E [0,7r/2]), 
and using that cos decreases on [0,7r], we conclude: 



6\ + ... + 0 n 

n 




This completes the proof of the Theorem. 

To prove the Corollary, we provide a standard argument, as done e.g. in ‘8], Theorem 33. 
Pick any initial state xq and consider the solution x(-) of the closed-loop system x = f(x,u — 
h(x)) with input 0 and x(0) = xo- Zero-reachability means that there is some finite-time input 
uq : [0, T] —> M such that, if zo(-) solves the closed-loop equations z = /(z, u — h{z)) with initial 
state Zo(0) = 0 and this input no on the interval [0,T], then zq(T) = x$. Consider now the 
input n obtained by the formula u(t) = uq(€) for t < T and u(t ) = 0 for t > T, and let z(- ) be 
the solution with initial state z(0) = 0 and this input n; by causality, z(t ) = ^o(t) for t < T, 
and hence z(T ) = xq = x(0), from which it follows that z{t + T) = x(t) for all t > 0. Showing 
x(t) —> 0 as t —> oo is the same as showing z(t) —> 0 as t —> oo. Let yi be the outputs of the 
subsystems when using input n (and zero initial state). Since u E L 2 and ||y|| < c||u|| < oo, 
we have that yi E L 2 for each of the intermediate outputs. Since each subsystem is input to 
state Testable, meaning that L 2 inputs (and zero initial state) produces L 2 state trajectories, 
we have that the complete state z is in L 2 . Finally, as z is a trajectory of a semiflow in finite 
dimensions, we must have that z(t) —* 0, by a Barbalat’s Lemma type of argument (see e.g. f2]). 

Finally, we review in the present context a weaker version that applies when = 2, basically 
part of the statement of the classical Passivity Theorem. Suppose that the first system is OSP 
but the second system is only known to be passive, in the sense that no estimate < 

72 ( 2 / 1 , 2 / 2 >t may hold, but we do know that ( 2 /i, 2 / 2 )r > 0 for all T > 0. Then, yo = — 2/2 implies 
that: 


II2/1II t ^ 7 i(« + 2 /o, 2 /i)t 

= 7 i(«, 2 /i)t ~ 71(2/2, 2/1 )t < 7 i(' u , 2 /i)t 

and so the system with output y\ is OSP, and in particular, L 2 stable. If, in addition, the second 
system is also L 2 stable, then stability to y 2 holds as well. 


5 Linear Systems 

The condition that a system be OSP is of course a restrictive one, but the concept of OSP 
system is thoroughly well-studied, and examples of passive systems abound, especially, but not 
only, for linear systems. We collect here some facts, mostly well-known, regarding the linear 
case. 

For a stable linear system with transfer function G(s), the secant gain can be characterized 
as the smallest 7 such that 

\G{iui)\ 2 < 7 ReG(*w) VcuEM. ( 6 ) 
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A proof is as follows. First of all, squaring the expression below and expanding (y — ( 7 / 2 ),?/ — 
( 7 / 2 ))ti one easily sees that the definition of OSP system is equivalent to the requirement that 

\\y - <Cl/2)u\\ T < ( 7 / 2 ) IMIt (7) 

for all i/o pairs and all T, which means 7 S is the smallest number such that the L 2 -induced 
norm of u 1 —» y — ( 7/2 )u is < 7 / 2 . For linear systems, induced L 2 -induced norm corresponds 
to H 0 o gain, that is to say, j s is the smallest number so that sup^gjj \G{iuj) — ( 7 / 2 )! < 7 / 2 . 
Writing \G(iuj ) — ( 7 / 2)| 2 = ( G{iu ) — ( 7/2 ))(G(ilu) — ( 7 / 2 )) and expanding, one has ©. 

An equivalent formulation of m is via the following analog of the estimate ©: 

\G{iu)\ < 7 cos 6{G(iu:)) Vw£R. ( 8 ) 

where we are denoting now by 9(p) the argument of a complex number p. Since G is analytic 
on Re A > 0 (stability), the maximum modulus principle for analytic functions implies that 
same estimate is obtained when maximizing not merely over A = iui purely imaginary, but also 
over all complex numbers with nonnegative real part. 

If we write G(s) = p{s)/q(s) as a quotient of two polynomials, condition © can be also 
written as 

\p(iu>)\ 2 < 7 Re \p{iuj)q(iuj)]. 

For example, for a one-dimensional system x = —ax + (3u with output y = x, the transfer 
function is /3/(s + a), so that p(ico) = /3 and Re \p{iuj)q{iui)\ = a/3 for any lo, from which it 
follows that 7 S = /3/a , and the classical result is obtained. On the other hand, as is well-known 
for OSP systems, G(s) must have relative degree at most one (the condition Re G(iu) > 0 is 
otherwise violated). Therefore, cascades, as studied here, of two or more such one-dimensional 
systems are not OSP themselves. 

For linear systems, a sufficient condition for a system to be OSP is that its transfer function 
G(s) be strictly positive real (SPR), meaning that G(s — e) is positive real for some e > 
0, or equivalently (see e.g. 0, Lemma 10.1) that it be stable (all poles have negative real 
part) and satisfy ReG(ia;) > 0 for all lo G R and lim^^oo cu 2 Re G{iuj) > 0. (Note that our 
transfer functions are strictly proper, by definition, since we are considering state-space systems 
with no direct i/o term; for non-strictly proper transfer functions, the condition is slightly 
different.) This provides a large class of examples; for instance, any transfer function of the 
form (s + a)/(s 2 + as + b) with b > 0 and 0 < a < 2 \fb is SPR if and only if 0 < a < a (0, 
Exercise 10.1). That SPR implies OSP can be proved using the Kalman-Yakubovich-Popov 
(KYP) Lemma. The converse implication does not hold: s/(s 2 + s + l) is not SPR, since it fails 
the test just quoted with (a = b = 1, a = 1/2) or just by noting that there is an imaginary axis 
zero, since ReG(O) = 0, but it is OSP, since \p(iuj)\ /Re [p(iuj)q{iui)] = 1 < 00. 

More generally, for not necessarily linear systems, if there exists some nonnegative definite 
smooth function V on states with the property that, for some 7 > 0, 

\7V(x).f(x,u) < —y 2 + 7 «!/ 

for all x £ M n , u £ M, and y = h(x ), then the system is OSP. Indeed, integrating along 
solutions corresponding to x(0) = 0, and using that V is nonnegative definite (so that R(0) = 0 
and V(x(T)) > 0), one has that 

0 < V(x(T)) - V'(O) 

[ T fT 

< — y(s) 2 ds + 7 / u(s)y(s)ds 

J 0 Jo 
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and thus || y||y < 7 (u,i/)t as claimed. This property can be checked by means of nonlinear 
versions of the KYP Lemma, see e.g. mm. 


Yet another way of stating the estimate © is in terms of integral quadratic constraints 
(IQC’s), cf. (S: one may equivalently write “ w T Mw > 0” in L 2 for i/o pairs w = ( u,y)' and 
where: 


M = 


( 0 7/2 \ 

V7/2 -l) 


The powerful tools for analysis of IQC’s, based on LMI’s, as developed by Megretski and Rantzer 
and others, should thus be useful for the study of secant gains. (We wish to thank R. Sepulchre 
for suggesting this reformulation.) 


We pointed out that the induced L 2 gain 700 is upper bounded by the secant gain j s . In 
general, one has the strict inequality 7^ < 7^. For example consider the linear system with 
transfer function 

2s+ 1 
s 2 + s + 1 ’ 


This is a scalar multiple of (s +1 /2)/(s 2 + s +1), so it is SPR by the criterion mentioned earlier, 
and hence OSP. Explicitly: 


and 


\p(iLo)\ 2 1 + 4 J 2 

sup--—= sup —- y 

weR Re \p{iuj)q(iu)] wgR 1 + w 


sup 

1 

= sup \ 

/ 1 + 4w 2 

1 — iui — u 2 

1 — UJ 2 + UJ 4 

ugM 

uGR ' 


= ^ 2 + (2/3)v^T ~ 2.25 < 4 


(the maximum value is achieved at ui = 1/2 \/ y/21 — 1). Graphically, we can see these conclu¬ 
sions from Figure [3 which shows that the smallest circle of the form |s — 7/2| < 7/2 which 



Figure 2: Nyquist plot for and circle |s — 2\ < 2, 7 S = 4, 7^ « 2.25 


contains the Nyquist plot must have 7 = 4 (circle shown), so that this is the value of 7 S , but 
7oo ~ 2.25 because the plot fits in a circle (not shown) centered at the origin with radius ~ 2.25. 

To conclude, let us provide a direct proof of the main theorem in the linear case. This 
proof, when specialized to linear one-dimensional systems, is basically the same as the proof 
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given in m- We assume a unity negative feedback about the cascade in Figure where each 
system has transfer function Gi(s ) and secant gain 7,. As remarked earlier, this means that an 
estimate |Gj(A)| < 7^ cos 0(Gj(A)) as in (JHJ) holds for every A with real part > 0. If the closed- 
loop were not to be stable, then there would exist a pole A with real part > 0 for G/( 1 + G), 
where G = G\ ... G n is the open-loop system. For any such A: 

Gi(A)G 2 (A)...G n (A) = -l (9) 

from which we conclude, writing 0^ = 0(Gj(A)), that 0i is a multiple of ir. Moreover, taking 
absolute values in © and using |Gj(A)| < 7jcos0j, we have also that 

1 < 7172 ... 7 n II” =1 cos 0i 

and so 71... 7 n > (sec^) n , again using the convexity of — Incos x. Thus no such poles can 
exist, if the hypothesis of the theorem holds. 
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